Disorder-mediated electron valley resonance in carbon nanotube quantum dots 
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We propose a scheme for coherent rotation of the valley isospin of a single electron confined 
in a carbon nanotube quantum dot. The scheme exploits the ubiquitous atomic disorder of the 
nanotube crystal lattice, which induces time-dependent valley mixing as the confined electron is 
pushed back and forth along the nanotube axis by an applied ac electric field. Using experimentally 
determined values for the disorder strength we estimate that valley Rabi oscillations with a period 
on the nanosecond timescale are feasible. The valley resonance effect can be detected in the electric 
current through a double quantum dot in the single-electron transport regime. 
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Introduction. The conduction and valence bands of 
graphene and carbon nanotubes (CNTs) form two valleys 
as the bands approach each other at two non-equivalent 
points (K and K') of the Brillouin zone pp. This twofold 
degeneracy of the electronic spectrum implies that the 
valley degree of freedom of an electron can be regarded 
as one bit of classical or quantum information, in analogy 
with the electron spin [2]. To evaluate the potential of 
encoding information in the valley degree of freedom, it 
is necessary to explore the physical mechanisms which 
could provide control over the valley state, as well as 
those leading to the loss of information. 

Recent theoretical proposals suggest the use of valley- 
polarized edge states of zigzag graphene nanoribbons in 
a valley filter device [2], and non-adiabatic magnetic 
field sweeps for manipulation of the valley state of a 
single electron in a graphene Aharonov-Bohm ring [3J. 
However, the uncontrolled and strong valley mixing due 
to edge irregularities of nanostructured graphene might 
pose a significant challenge towards the realization of 
these ideas. This difficulty motivates the study of valley 
physics in CNTs, which have a rolled-up and therefore 
edge- free geometry. In fact, high-quality CNT quantum 
dots (QDs) have been fabricated recently [3HS] and the 
fourfold (spin and valley) quasi-degeneracy of the QD en- 
ergy levels as well as very weak disorder have been exper- 
imentally confirmed [4j [5] . These ultraclean QD devices 
offer a chance for experimental realization of spin control 
via spin-orbit interaction [7] and manipulation of a com- 
bined spin- valley qubit (the "Kramers qubit" ) utilizing a 
magnetic field, valley mixing, and bends in a CNT [8]. 

In this work we propose a scheme for coherent control 
of the valley isospin of a single electron, which is confined 
in a QD in a straight CNT at zero magnetic field. The 
effect, which we name electron valley resonance (EVR) in 
analogy to electron spin resonance (ESR) , relies on short- 
range atomic disorder (substitutionals, adatoms, vacan- 
cies) of the CNT crystal lattice. Disorder induces har- 
monic time-dependent valley mixing as the confined elec- 
tron is pushed back and forth along the CNT by an ap- 




FIG. 1: (Color online) Schematic of the valley resonance 
setup, (a) Two gate electrodes close to the CNT confine an 
electron to the QD region, (b) The confinement potential 
V(z) is periodically modulated in time by an ac voltage on 
the left gate, (c) The modulation of the confinement potential 
is followed by a displacement of the electron density n(z) . The 
interplay of this driven motion and the disorder on the CNT 
induces ac valley mixing leading to the valley resonance effect. 



plied oscillating (ac) electric field (Fig. [I]), which in turn 
induces Rabi oscillations between the two valley states. 
The ac electric field can be induced by applying an ac 
voltage component on one of the confining electrodes, as 
demonstrated in single-spin control experiments in GaAs 

qds ma. 

Using a microscopic model [HJ [12] of disorder-induced 
valley mixing, we map the valley dynamics of the oscil- 
lating electron to the two-level ESR problem. Together 
with a statistical treatment of possible disorder config- 
urations, we find that the ratio between typical energy 
scales of the ac and static valley mixing is Z/L, where 
Z is the displacement amplitude of the electron and 2L 
is the width of the electronic wave function [see Eq. Q]. 
As the static valley- mixing energy is measurable [HlS]> 
our finding can be used to estimate the ac component 
and hence the frequency of the valley Rabi oscillations. 
For a numerical example (see below) using data from [4] , 
with Z = 2 nm and L = 50 nm we find that the time 
needed for a half Rabi oscillation between the two valley 
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states is ~ 1.6 ns. We also describe a setup that could be 
used to detect the EVR effect via a measurement of the 
electric current in the regime of single-electron transport 
through a CNT double QD (DQD). 

Microscopic theory of EVR. To model the electronic 
states we use the envelope function approximation and 
take the coordinates defined in Fig. [I] The Hamilto- 
nian includes four terms: (i) the kinetic energy [TTj 
-^kin = vpfoaiPx + o~2Pz) with Fermi velocity v-p, cir- 
cumferential (longitudinal) electron momentum (Pz)i 
and o~i (tj) being Pauli matrices acting in sublattice (val- 
ley) space; (ii) curvature-enhanced spin-orbit interaction 
[HHSHIS] -ffso = t 3 s z (Ao + ctiAi) with axial spin com- 
ponent s z and chirality-dependent on-site (off-site) spin- 
orbit matrix element A (Ai); (iii) QD confinement po- 
tential Vq(z), and (iv) disorder induced by atomic defects 
[11) which randomize the on-site energies on the crystal 
lattice, 



H dis (r) = fl ccll ^2 ^2 cr A T o + T r {'fii < j)}Ui C rS(r-ri lT ), 

l <j=A,B 

(1) 

where Ui a is the on-site energy on site a g (A, B) = 
{+,—) in unit cell I. Here Sl co ii is the area of the unit 
cell of the graphene lattice, o~a,b = (co ± o"3)/2, T r (ip) = 
cosipTx + sin</jr 2 , (pi a = a (2K ■ r la - -q) + ^ B 27r/3, r\ 
is the chiral angle of the CNT, r\ a is the position of 
the lattice site la, and K is the vector pointing to the 
K point of the Brillouin zone. The presence of the off- 
diagonal valley operator r r in H<n s reflects the fact that 
atomic disorder, due to its short-range character, allows 
for large momentum transfer upon scattering, including 
intervalley (K ■<-» K') transitions. Our valley rotation 
scheme relies on this K-K' coupling. 

EVR is induced by pushing the confined electron back 
and forth along the z axis by applying an ac voltage 
on one of the gates (Fig. [lj. We model this time- 
dependent displacement with the confinement poten- 
tial V(z,t) = Vo(z — Z smut) with the ac frequency 
uj/2it. To map the system to a two- level ESR problem, 
we first transform the complete Hamiltonian into the 
reference frame co-moving with the displaced electron. 
Thus we perform the time-dependent unitary transfor- 



mation U(t) 



p iZ sin(ujt) 



P*l h , which leaves the position- 



independent kinetic energy and spin-orbit Hamiltonians 
i/kin and H so invariant, but renders the confinement po- 
tential time-independent V(z,t) t— > Vq(z) and the disor- 
der term time-dependent -ffdis(' r ) ^ H dis (x, z + Z sin uot) . 
The additional term iW^JJ- 1 ^), arising from the time- 
dependence of the transformation, can be neglected as it 
does not induce transitions within the ground state QD 
level. 

As the next step, we restrict our consideration to the 
fourfold-degenerate (spin and valley) QD ground-state. 
For our purposes we can estimate the corresponding 
eigenfunctions of Hy i i n + Vo(z) as Gaussians, having equal 



weight on the two sublattices: ^k,s(x, z) = \K)x s — 
^(Gz(z),G z (z),0,0)x,, and * K ,, a (x,z) = \K') X . = 
^=^(0,0,Gl(z),G l (z))xsi with spinors \s having spin 
projection s e (1\4-) = (+> ~~ ) along the z axis, cir- 
cumferential wave number q > 0, CNT radius R, and 

-z 2 /2L 2 

Gl(z) — ^.1/4^2^ ' Projecting the spin-orbit Hamilto- 
nian to this four-fold degenerate subspace yields H so = 
lA so s z {\K')(K'\ - \K){K\), where A so = 2(A + Aj); 
experimentally reported values [H [T5] of Aso are in 
the range 0.17-2.5 meV. The form of H so implies that 
spin-orbit interaction induces an energy splitting Ago 
between two valley states having the same spin direc- 
tion. Therefore, as in the case of ESR in CNT QDs 
7J, no magnetic field is needed for EVR . Henceforth we 
restrict our considerations to the limit of small displace- 
ments Z/L <C 1. By projecting the disorder Hamiltonian 
to the four-dimensional subspace of interest we find 



Ha 



(be**- 



iut)\K')(K\+h.c, (2) 



where valley-diagonal terms are omitted. The real quan- 
tites b, <f>, 6 ac and </> ac describe disorder-induced static 
and ac valley mixing and can be expressed in terms of 
Uicr and fy VtS - Randomness of the disorder configuration 
Ui<j implies the randomness of those quantities as well. 
Assuming a homogeneous and uncorrelated distribution 
of the atomic defects ^UioUva 1 ) = {Uf^&Wa't w ith zero 
average (£/j CT ) = 0, we find (b) — (b ac ) = and 
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Equations ^ and Q are the central results of this 
work. The ac valley-mixing term oc 6 ac in Eq. ([2]), in- 
duced by the simultaneous presence of the ac electric 
field and atomic disorder, allows for coherent rotations 
of the valley isospin similarly to a transverse magnetic 
field in ESR [21]. Furthermore, using the remarkably 
simple expression in Eq. Q we can estimate the cor- 
responding Rabi frequency from a measurement of the 
static valley-mixing matrix element b. For example, in 
@] a spin-orbit splitting of Ago = 370 ^eV and a val- 
ley gap of Akk ' = 65 /zeV was found. Identifying Akk> 
with 2yJ (b 2 ), and taking L — 50 nm and Z — 2 nm, from 
Eq. Q we find an estimate for the strength of the ac 
valley-mixing term in this sample as U (b 2 c ) = 1.3 fieV. 
This value translates to a ir 

h /V( b lc) ~ 1-6 ns long half 
Rabi cycle between two orthogonal valley states at reso- 
nant driving fno = \/Al Q + Ab 2 . 

Detection. We now describe a setup where EVR could 
be detected (Fig. [2]). A serially coupled DQD between a 
source and a drain lead is tuned appropriately (see below 
for details) , such that during the single-electron transport 
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FIG. 2: (Color online) Schematic of EVR detection. The 
disorder-mediated ac valley mixing drives transitions between 
the states on the left QD. Interdot tunneling between K' L and 
K R has the amplitude re"* T = t(K R \K' L ). 



process the electron can be trapped in a specific valley 
state in the left QD. (\K L ,0) in Fig. [2). The valley 
state of a trapped electron is changed due to the EVR 
mechanism when an ac electric field is applied to the 
left QD. This allows the particle to exit via the ground 
state level of the right QD, thus EVR is detected via the 
measurement of the electric current through the DQD. 
Since this method is based on single-electron transport, 
it is unaffected by strong correlations, unlike the Pauli 
blockade effect E2 El QU . 

To describe the transport process outlined above, 
from now on we consider only spin-f electrons; 
hence the spin-orbit Hamiltonian simplifies to _ff so = 
iA so (\K')(K'\ - \K){K\). The description of the spin- 
1 electrons is completely analogous. As shown in Fig. 
[2j the energy levels in the DQD are tuned by gate elec- 
trodes so that the higher-lying sublevel in the left QD 
is aligned with the lower-lying sublevel in the right QD. 
We consider single-electron transport via the (0, 0) — > 
(1,0) — > (0,1) — > (0,0) transport cycle where (n,m) 
refers to the charge configuration with n (to) electrons in 
the left (right) QD. Interdot tunneling is assumed to be 
spin- and valley-conserving, therefore the static Hamil- 
tonian of a single electron in the DQD is 
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where we use the basis \K,0), \K',0), \0,K), \0,K'), and 
L, R = 0, K, K' in \L, R) refer to the occupation and val- 
ley state in the left and right QDs. The left-right energy 
detuning, needed for the the level alignment described 



above, is e = (a/A 
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For the 



sake of simplicity and clarity, from now on we assume 
bL,B, <C Aso (as found in recent experiments [HE]), and 
t <C Aso (tunability of the interdot tunneling has been 
demonstrated in CNT DQDs [19]). A crucial point for 
our detection scheme is that the electron interacts with 
a different set of impurities in the two different QDs, and 
therefore 6l and 4>l are in general different from b R and 



4>r. Together with the condition b^ : R <C Aso, this im- 
plies that the higher-lying energy eigenstate in the left 
QD \K' L ) has a small K component, the lower-lying state 
in the right QD \K R ) has a small K' component, and 
\K' L ) is typically not orthogonal to \K R ). This leads to 
a finite tunneling matrix element Te % ^ T — t(K R \K' L ) be- 
tween the two aligned levels \K' L ,0) and \0,Kr). This 
nonzero matrix element is required for the EVR detection 
scheme we outline below. 

The level structure shown in Fig. [2]results in a blockade 
effect in a finite source-drain voltage bias. If the incoming 
electron enters the QD in the higher-lying | K' L , 0) state 
then it can move through the DQD easily because of the 
complete hybridization of \K' L ,0) with \0,K R ). How- 
ever, the condition t <g; Aso ensures that hybridization 
of the lower-lying state \Kl,0) with (0,1) states is small 
(< tj Aso m amplitude); therefore an electron occupying 
that state blocks transport as it exits the DQD slowly. 
Assuming equal tunneling rates T at source and drain, 
from a transport model based on the secular approxi- 
mation (T <C t, bL,R, Aso) of the Born-Markov master 
equation we estimate that the current through the DQD 
without ac driving is much smaller than barrier trans- 
parency r, 



I /e ps 2r(VA so ) 2 « r. 
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This transport blockade allows for detection of the 
EVR effect described above. If an ac electric field is ac- 
tive in the left QD and pushes back and forth the electron 
within the left dot, occupying the blocking lower-lying 
state 0) in Fig.[2j then it undergoes a Rabi transition 
to the higher- lying \K' L ,0) state, which allows it to exit 
the DQD to the drain due to the strong hybridization of 
\K' L , 0) and |0, Kr). This process is most effective around 
resonant driving, i.e., when hui w \J Ag Q + 46^. In fact, 
below we show that around the resonance condition and 
an appropriate tuning of the tunneling amplitude r the 
current through the DQD approaches its maximal value 
2er/7. 

To provide an analytical result for the ac-field-induced 
current we neglect perturbative hybridization amplitudes 
between (1,0) and (0,1) states which are < i/Aso- In 
this simplified picture the higher-lying level 1 0, K' R ) on 
the right QD is not involved in the transport problem as 
it is not hybridized with (1,0) states. After a transfor- 
mation into the rotating frame in the left QD and the 
counter-rotating frame in the right QD, and within the 
rotating wave approximation (RWA), the dynamics in the 
remaining three-level system can be described by 
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with the ac valley-mixing term oc 6 ac and the detuning 
parameter 5 — hco — \J A| Q + Ab 2 L . The phases 4> ac and 
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FIG. 3: (Color online) Current as function of the dimension- 
less detuning 8 = 5/6 ac of the driving frequency from reso- 
nance, for different values of the dimensionless tunnel ampli- 
tude f = r/6 ac . Inset: Energy diagrams corresponding to the 
two indicated points of the main plot. 



4> T have been eliminated by a specific choice of the basis 
states. A small term oc 6 ac 6/Ago has been neglected in 
Eq. ([7]). The leading-order (in V) analytical result for the 
current, obtained from the Born-Markov master equation 
transport model is 
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Using Eq. ^ , in Fig. [3] we plot the current as a func- 
tion of dimensionless detuning 6 = 5/b ac for three dif- 
ferent values of the dimensionless tunneling amplitude 
f = r/6 ac . The three curves correspond to three dif- 
ferent regimes: (i) For f = 6.0 the current shows two 
Lorentzian peaks as the function of detuning, which have 
width 6 ac /\/2 and maximal current er/4 in the entire 
f >1 regime. In this situation, hybridization of \K' L ,0) 
and \0,Kr) results in a 'bonding' and an 'antibonding' 
state having energies =Fr, depicted as dashed lines in the 
left diagram of the inset of Fig. [3j The ac field can be 
resonant with only one of these states (in the inset it is 
the bonding one) , as the current peak width is set by the 
Rabi energy 6 ac , which is significantly smaller than the 
peak separation 2r in this case. This parameter regime 
seems to be suitable for EVR detection as the resonant 
current is large, i.e., comparable to the tunneling rate 
of the lead-dot barriers (see, however, estimate of the 
background current Jo below), (ii) For f = 0.7, the en- 
ergy separation 2r between the two peaks of regime (i) 
becomes comparable to the Rabi energy 6 ac , and there- 
fore the two peaks merge into a single one. This pa- 
rameter setting is optimal for EVR detection as the peak 
value of the current is 2er/7, exceeding that in regime (i). 
(iii) For f = 0.2, the Rabi energy 6 ac becomes dominant 
over interdot tunneling r. This implies that the mini- 
mal energy distance between the two (1,0) states cannot 
drop below b ac , and current is maximal if they hybridize 
equally (oc r/6 ac ) with |0, Kp) (this maximal-current sit- 



uation is shown in the right diagram of the inset of Fig. 
[3). This f <C 1 regime is less suitable for EVR detection 
as the peak current is only a small fraction 2f 2 T of the 
barrier transparency V . The line shape is approximately 
a Lorentzian for small detuning, with width & a e/v5- 

To take a numerical example, we use again b ac — 
1.3 fieV, A so = 370 /ieV, b L = b R = 32.5 ^eV, cj> R -cj> L = 
7r/3, which yields r s» 0.075t. This implies that by tun- 
ing t to 90.2 [ieV, 10.5 /xeV and 3.0 /zeV, one arrives to 
the previously plotted and described situations (i), (ii) 
and (iii), respectively. Furthermore, it is reasonable to 
assume a lead-dot tunneling rate T = 62.5 MHz, which 
translates to a measurable current scale er 10 pA, and 
a tunneling-induced level broadening of hT w 40 neV 
<C b ac ,r in all cases. We also estimate the EVR detec- 
tion visibility by comparing the background current Jo 
[see Eq. dfj|], i.e., the current in the absence of the ac 
electric field, to the current at resonant EVR driving. 
In the cases (i), (ii) and (iii), the background current Jo 
is estimated from Eq. ^ as 1.2 pA, 16 f A and 1 fA, 
respectively. For this numerical example, we can con- 
clude that in cases (ii) and (iii) the background current 
is much smaller than the EVR contribution due to the ac 
electric field (~ 1 pA from Fig. [3] and er = 10 pA), and 
therefore EVR could indeed be detected in this transport 
measurement. 

Valley relaxation, e.g., mediated by valley mixing and 
accompanied by phonon emission, can affect the dynam- 
ics of the EVR. By including valley relaxation in our 
transport model, we have found that the results shown 
in Fig. [3] do not change significantly as long as the valley 
relaxation rate is below b ac /h and r/h. In principle, EVR 
could also be detected by generalizing the proposed setup 
to pulsed-gated transport or charge sensing experiments 
[SI OH [ini HO] , which may allow for measuring the valley 
decoherence time scale. Finally, we note that the ideas 
described in this paper in the context of CNTs could be 
generalized to QDs in other multivalley materials, such 
as silicon. 
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